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Abstract 

We calculate dark matter scattering rates in the minimal supersymmetric extension 
of the Standard Model (MSSM), allowing the soft supersymmetry-breaking masses of the 
Higgs multiplets, m 12 , to be non-universal (NUHM). Compared with the constrained MSSM 
(CMSSM) in which m 12 are required to be equal to the soft supersymmetry-breaking masses 
mo of the squark and slepton masses, we find that the elastic scattering cross sections may 
be up to two orders of magnitude larger than values in the CMSSM for similar LSP masses. 
We find the following preferred ranges for the spin-independent cross section: 10~ 6 pb 
^ csi ~ 10~ 10 pb, and for the spin-dependent cross section: 10~ 3 pb ^ o~sd, with the 
lower bound on asi dependent on using the putative constraint from the muon anomalous 
magnetic moment. We stress the importance of incorporating accelerator and dark matter 
constraints in restricting the NUHM parameter space, and also of requiring that no unde- 
sirable vacuum appear below the GUT scale. In particular, values of the spin-independent 
cross section another order of magnitude larger would appear to be allowed, for small tan/3, 
if the GUT vacuum stability requirement were relaxed, and much lower cross-section values 
would be permitted if the muon anomalous magnetic moment constraint were dropped. 
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1 Introduction 



There have been many previous studies of the elastic scattering rates of supersymmetric 
relic particles on matter in the minimal supersymmetric extension of the Standard Model 
(MSSM) [1] - [6], assuming conservation of R = (— i) 3B + L + 2S ' j where B is the baryon number, 
L the lepton number and S the spin, so that the lightest supersymmetric particle (LSP) is 
absolutely stable. As in most previous studies, we assume this to be the lightest neutralino 
X [7]. In this paper, we refine and extend previous calculations of the elastic scattering 
rates when the input soft supersymmetry-breaking scalar masses for the Higgs multiplets 
are allowed to be non-universal at the input GUT scale, the non-universal Higgs model 
(NUHM). 

As we discuss later in more detail, it is important to impose the constraints due to 
accelerator experiments, including searches at LEP, b — > 57 and (optionally) the muon 
anomalous magnetic moment, — We assume also that most of the cold dark matter is 
composed of LSPs, with relic density 0.1 < Q x h 2 < 0.3, while being aware that the lower 
part of this range currently appears the most plausible [8]. 

In the constrained MSSM (CMSSM), in which all the soft supersymmetry-breaking scalar 
masses m are assumed to be universal, including those for the Higgs doublets #1,2, the 
underlying parameters may be taken as m , the soft supersymmetry-breaking gaugino mass 
TO1/2 which is assumed to be universal, the trilinear supersymmetry-breaking parameters A 
that we set to zero at the GUT scale in this paper, and the ratio tan/5 of Higgs vacuum 
expectation values. In the NUHM [9] - [13], there are two additional free parameters, the two 
soft Higgs masses or equivalently the Higgs superpotential coupling /x and the pseudoscalar 
Higgs boson mass uia- These would be fixed by the electroweak symmetry-breaking vacuum 
conditions in the CMSSM, up to a sign ambiguity in /i, in terms of the other parameters 
(mo, m\/2i A , tan/3). We use the parameters (mo, my 2, m A, Aq, tan/3) to parametrize the 
more general NUHM. 

As we have pointed out previously [12, 13], this six-dimensional NUHM parameter space is 

significantly restricted by the requirement that no undesirable vacuum appears when one uses 

the renormalizat ion-group equations to run the soft supersymmetry-breaking parameters 

between the input GUT scale and the electroweak scale. This requirement constrains the 

non-universalities of the Higgs masses: m 8 = sign(mf)|mj/mo| : % — 1,2, which in turn 

restricts the range of elastic scattering cross sections that we find 1 . 

1 This extended stability requirement would also exclude non-universalities for the input squark and 
slcpton masses that allowed their physical values to be similar, an assumption that might lead to much 
larger elastic scattering cross sections than we find here. 
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The allowed regions in the (fj,,m A ), (//, M 2 ) and (mi/2,mo) planes for certain discrete 
values of the other NUHM parameters have been described in [13]. Our first step in this paper 
is to provide contours of the elastic scattering cross sections in selected planes, providing a 
first comparison with the CMSSM points that appear in these planes. Secondly, we display 
the ranges of elastic scattering cross sections that are allowed in these planes, as functions of 
the LSP mass. In general, we find that the spin-independent elastic scattering cross sections 
may be up to two orders of magnitude larger than values in the CMSSM for similar LSP 
masses, and another order of magnitude larger if the GUT vacuum stability requirement is 
relaxed. Thirdly, we display ranges of the elastic scattering cross sections as functions of the 
LSP mass for all allowed values of the other NUHM parameters. 

We review the NUHM in Section 2, including the experimental and phenomeno logical 
constraints on its parameter space. Then, in Section 3, we summarize our treatment of the 
elastic scattering matrix elements and display contours of the cross sections in various planar 
projections of the NUHM parameter space. Section 4 presents and discusses the ranges of 
the cross sections attainable in the NUHM. Finally, Section 5 draws some conclusions from 
our analysis. 

2 The NUHM and Constraints on its Parameter Space 

We assume that the soft supersymmetry-breaking parameters are specified at some large 
input scale M x , such as the supergravity or grand unification scale. Motivated by restrictions 
on flavour-changing neutral interactions, we assume that squarks and sleptons with the 
same Standard Model quantum numbers have universal soft supersymmetry-breaking scalar 
masses at this input scale. With the weaker justification provided by some GUTs, we further 
assume universality between the soft scalar masses of squarks and sleptons. However, in 
the NUHM studied here, we allow the soft supersymmetry-breaking scalar contributions 
to the masses of the Higgs supermultiplets at M x to be free non-universal parameters. 
Their running from M x down to low energies relates m\(M x ) and m\(M x ) to the Higgs 
supermultiplet mixing parameter /i and the pseudoscalar Higgs mass uia- Therefore, we use 
as free parameters /i(m z ) = ji and itia(Q) = tua, where Q = (m^m^) 1 / 2 , in addition to the 
parameters (m (M x ), m 1/2 (M x ), A , tan (3) used in the CMSSM 2 . 
The electroweak vacuum conditions may be written in the form: 

ml(Q) = mj(Q) + m 2 (Q) + 2/i 2 (Q) + A A (Q) (1) 

2 In this paper, we use the input A a = for definitcncss, noting that the range of effective low-energy 
values of A, after renormalization below the GUT scale, is quite limited. 
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and 

2 m\ -m 2 tan 2 /3 + |m|(l -tan 2 /3) + 
" = tan^-l + A^ ' (2) 

where A A and A^ 1 ' 2 ) are loop corrections [14, 15, 16] and mi j2 = rni^{m z ). We incorporate 
the known radiative corrections [14, 17, 18] ci,c 2 and c M relating the values of the NUHM 
parameters at Q to their values at m z : 

m\(Q) = m\ + ci 

m\{Q) = ml + c 2 

» 2 (Q) = /i 2 + c M . (3) 

Solving for m\ and m|, one has 

m 2 (l + tan 2 /3) = m\{Q) tan 2 /3 - /i 2 (tan 2 /3 + 1 - A^) - (c x + c 2 + 2c M ) tan 2 /3 

-A A (Q) tan 2 /3 - ^m 2 (1 - tan 2 (3) - A« (4) 

and 

m 2 (l + tan 2 /5) = m 2 A (Q) — /i 2 (tan 2 /3 + 1 + A^) — (ci + c 2 + 2c M ) 

-A A (Q) + im| (1 - tan 2 /?) + A«, (5) 

which we use to perform our numerical calculations. 

It can be seen from (4) and (5) that, if m A is too small or /i is too large, then m\ and/or 
m\ can become negative and large. This could lead to m\(M x ) + n 2 (M x ) < and/or 
m\{Mx) + n 2 (M x ) < 0, thus triggering electroweak symmetry breaking at the GUT scale. 
The requirement that electroweak symmetry breaking occurs far below the GUT scale forces 
us to impose the conditions m\{M x ) + fi(M x ),ml(M x ) + fi(M x ) > as extra constraints, 
which we call the GUT stability constraint 3 . We have discussed in [13] issues related to the 
NUHM renormalization group equations (RGE's) and their solutions. 

We assume that R parity is conserved, so that the LSP is stable and is present in the 
Universe today as a relic from the Big Bang. Searches for anomalous heavy isotopes tell 
us that the dark matter should be weakly-interacting and neutral, and therefore eliminate 
all but the neutralino and the sneutrinos as possible LSPs. LEP and direct dark-matter 
searches together exclude a sneutrino LSP [20], at least if the majority of the CDM is the 
LSP. Thus we require in our analysis that the lightest neutralino be the LSP. We include in 
3 For a different point of view, however, see [19]. 
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our analysis all relevant coannihilation processes involving this LSP and sparticles that may 
become near-degenerate in various regions of the NUHM parameter space. We restrict our 
attention to regions of the NUHM parameter space where 0.1 < Vl x h 2 < 0.3. 

We impose in our analysis the constraints provided by direct sparticle searches at LEP, 
including that on the lightest chargino x ±: m x ± ~ 103.5 GeV [21], and that on the selectron 
e: m~ e <: 99 GeV [22]. Another important constraint is provided by the LEP lower limit on 
the Higgs mass: m H > 114.4 GeV [23] in the Standard Model 4 . The lightest Higgs boson h 
in the general MSSM must obey a similar limit, which may in principle be relaxed for larger 
tan/3. However, as we discussed in our previous analysis of the NUHM [12], the relaxation in 
the LEP limit is not relevant in the regions of MSSM parameter space of interest to us. We 
recall that rrih is sensitive to sparticle masses, particularly mi, via loop corrections [24, 25], 
implying that the LEP Higgs limit constrains the NUHM parameters. We also impose the 
constraint imposed by measurements of b — > S7 [26], as discussed in [13]. 

We take an agnostic attitude towards the latest value of the anomalous magnetic moment 
of the muon reported [27] by the BNL E821 experiment. The world average of a M = — 
now deviates by (33.7 ± 11.2) x 10~ 10 from the Standard Model calculation of [28] using e + e~ 
data, and by (9.4±10.5) x 10~ 10 from the Standard Model calculation of [28] based on r decay 
data. On some of the subsequent plots, we display the formal 2-cr range 11.3 x 10~ 10 < Sa^ < 
56.1 x 10~ 10 . However, in view of the chequered history of the theoretical Standard Model 
calculations of a M , we do not impose this as an absolute constraint on the supersymmetric 
parameter space. 

The results of applying the above constraints to various two-dimensional projections of 
the NUHM parameter space were described in [13]. 

3 Contours of the Cross Sections for Elastic Scattering 

The code we use to calculate the spin-independent and -dependent elastic dark matter scat- 
tering cross sections <Jsi,sd was documented in [2, 3], together with the ranges of values 
of the hadronic matrix elements that we use. The cross sections for protons and neutrons 
are similar within the quoted uncertainties in these matrix elements. There are other codes 
available [29] that include additional contributions to the scattering matrix elements, but a 
comparison [4] showed that the improvements were not essential for the CMSSM, and we 
believe they may also be neglected for our comparisons of the NUHM. 

4 In view of the theoretical uncertainty in calculating m/,, we apply this bound with just three significant 
digits, i.e., our figures use the constraint m,h > 114 GeV. 
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In [13], we analyzed NUHM dark matter in two ways: (i) fixing tan/3 = 10 and /i > 0, 
but choosing different values of \x and tua, rather than assuming the CMSSM values, and 
(ii) varying tan/3 for representative fixed values of \i and tua- We presented in [13] three 
types of slices through the NUHM parameter space, along (mi/ 2 ,m ) planes, (^m^) planes 
and (fi, M 2 ) planes. 

In this paper, we concentrate first on a few specific examples of these slices, presenting 
later more general results. We choose two representative examples each of the (TOi/ 2 ,m ) 
planes, (//, m^) planes and (//, M 2 ) planes shown previously. As we discuss later in more 
detail, the dependences of the cross sections on tan (5 are weaker than those on some other 
parameters, so we concentrate on planes for tan/5 = 10. The examples we choose are the 
(m 1/2 ,m ) planes for /j, = 400 GeV, m A = 400 GeV and /i = 700 GeV, m A = 700 GeV, 
corresponding to Figs. 2(a) and (d) of [13], the (//, m^) planes for m = 100 GeV, toi/ 2 = 
300 GeV and m = 300 GeV, mi/ 2 = 300 GeV, corresponding to Figs. 4(a) and (c) of [13], and 
the (/!, M 2 ) planes for m = 100 GeV, m A = 300 GeV and m = 300 GeV, m A = 500 GeV, 
corresponding to Figs. 8(a) and (c) of [13]. 

3.1 Examples of (mi/2,mo) Planes 

We display in Fig. 1 contours of (a,b) the spin-independent and (c,d) the spin-dependent 
elastic scattering cross sections, in the cases tan/3 = 10, and (a,c) /i = 400 GeV and m A = 
400 GeV, and (b,d) \x = 700 GeV and m A = 700 GeV. We assume here and in the subsequent 
figures that A — 0, m t — 175 GeV and m^m^)^ = 4.25 GeV. Here and elsewhere, the 
thickest contours denote decades in the cross-section values in pb, labelled by their exponents. 
The medium and thinnest lines are intermediate contours in the cross-section values, namely 
2x and 5x decades, as labelled. 

We notice immediately that the cross-section contours are nearly vertical at large m^, 
reflecting the fact that they become almost independent of m in the NUHM. We also notice 
that, within the GUT stability range (inside the black dot-dashed curves), the cross sections 
increase with m^. This is because the LSP becomes more Higgsino-like as mi/ 2 increases. 
However, the cross sections do decrease again for very large mi/ 2 beyond the GUT stability 
limit as the low energy scalar masses increase with mi/ 2 as does the light Higgs mass (though 
slowly) . 

The GUT stability requirement imposes mi/ 2 ^ 600(1100) GeV in panels (a,c) and (b,d), 
respectively. Because of the increases in the cross sections with mi/ 2 , there are in turn upper 
bounds on the cross sections, that would not be respected if GUT stability were disregarded. 
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tan (3 = 10 , p. = 400 GeV, m A = 400 GeV tan (3 = 10 , jx = 700 GeV, m A = 700 GeV 




m 1/2 (GeV) m 1/2 (GeV) 



Figure 1: Contours of (a,b) the spin-independent and (c,d) the spin-dependent cross sections 
(solid black lines) are superimposed on the NUHM (m^, mo) planes for tan/3 = 10 and 
(a,c) /i = 400 GeV and m A = 400 GeV, and (b,d) fi = 700 GeV and m A = 700 GeV. The 
near-vertical (red) dot-dashed lines are the contours m/, = 114 GeV, and the near-vertical 
(black) dashed lines at lower mi/ 2 are the contours m x ± = 103.5 GeV. The dark (brick red) 
shaded regions is excluded because a charged particle is lighter than the neutralino, and the 
darker (dark blue) shaded regions is excluded because the LSP is a sneutrino. The light 
(turquoise) shaded areas are the cosmologically preferred regions with 0.1 < Q x h 2 < 0.3. The 
dark (black) dot-dashed lines indicate the GUT stability constraint. There are two such lines 
for each panel and only the areas in between are allowed by this constraint. 
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Because the — 2 constraint would provide even stronger upper bounds on TO1/2, it would 
also impose stronger upper bounds on the cross sections. In Figs. 1 - 3, we have left off the 
<7 M — 2 contours to avoid confusion with the cross section contours we are highlighting here. 
For the case with \x = 400 GeV and m A = 400 GeV, the — 2 constraint places an upper 
limit on mi/ 2 of about 400 GeV and for the case with fx = 700 GeV and m A = 700 GeV, the 
limit is mi/2 ^ 450GeV. 

As already remarked, the cross sections themselves do not vary greatly with tan (3, but 
the interplay of the other constraints is rather tan/3-dependent. In particular, at large tan/5 
the g^ — 2 constraint would not reduce significantly the upper bounds on the cross sections. 

3.2 Examples of (//, tjia) Planes 

We display in Fig. 2 contours of (a,b) the spin-independent and (c,d) the spin-dependent 
elastic scattering cross sections, in the cases tan/3 = 10, and (a,c) mo = 100 GeV and 
— 300 GeV, (b,d) itlq = 300 GeV and mi/2 = 300 GeV. We see that there are large 
suppressions in the spin-independent cross section for \x ~ —100 GeV and m A ^ 500 GeV, 
reflecting a cancellation in the matrix element. Apart from this, the cross sections generally 
decrease with increasing and (to a lesser extent) tjia- In this sense, the lower bounds on 
I and rriA set an upper bound on the cross section in the allowed region. 
In the CMSSM, cancellations which drive the cross section to extremely small values 
occur only in the spin-independent case, and only for \x < 0. In the NUHM model, however, 
we find that there is a new source for a cancellation which affects the spin-dependent cross 
section for either sign of //. The reason this occurs is as follows. In the CMSSM, the 
spin-dependent cross section is dominated by the up-squark exchange term and, despite the 
difference in the relative signs of the up-type and down-type contributions, the total cross 
section remains reasonably large. On the other hand, in the NUHM case considered here, 
low-energy sfermion masses are affected by the splitting between the soft supersymmetry- 
breaking Higgs masses at the GUT scale, S = gfimf — ml) /4. When S 7^ 0, the up squarks 
get somewhat heavier and the down squarks somewhat lighter as S increases (see [13] for 
details on the effects of S in the renormalizat ion-group equations), opening up the possibility 
of a cancellation between the contributions. In Fig. 2(c,d), the cancellation occurs between 
the two 10~ 7 pb contours very close to, but outside the GUT stability curve. The GUT 
stability requirement, that bounds \fi\ ^ 700 GeV in the figures displayed, therefore provides 
lower bounds on the spin-dependent cross sections. These are somewhat lower than the values 
found at the CMSSM points indicated by crosses in Fig. 2. The same cannot be said for the 
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spin-independent cross-sections, because the very small cross sections due to cancellations 
occur within the GUT stability region, except when one applies the — 2 constraint, which 
excludes the /x < region. The interplay of the other constraints is more complicated: as 
usual, /i < is disfavoured by the m h , b — > 57 and — 2 constraints. The upshot for /i > 
is that the cross sections are bounded above by the & x h 2 constraint, so that they cannot be 
much more than a factor of ~ 10 greater than the CMSSM values. 

For the values of parameters chosen in Fig. 2 (a,c), we obtain a spin-independent cross- 
section of a S i — 2.6 x 1CT 9 pb and a spin-dependent cross-section of a SD = 5.4 x 1CT 6 pb in 
the CMSSM for /j, > (cf. the position of the crosses in the figures). In the NUHM, we find 
that the range of possible cross sections (for this case) is 3 x lCT 10 pb < <j S i ^3x 10~ 8 pb 
and 3 x 10~ 8 pb < a SD < 1.6 x 10~ 4 pb when all constraints other than — 2 are included. 
For the parameters in Fig. 2 (b,d), we find the CMSSM spin-independent cross section a si = 
1.9 x 1CT 9 pb, while the spin-dependent cross section is relatively unchanged. Note, however, 
that this CMSSM point would be excluded due to an excessive value for Vth 2 ( > 1). In the 
NUHM, this parameter choice is allowed and gives the range lCT 9 pb < a SI ^ 8 x 10~ 8 pb 
and 9 x lCT 9 pb < a S D ^ 2x lCT 4 pb for the elastic cross sections. In this case, there is 
essentially no 'bulk' cosmological region, and the spread in a is due to the region where the 
relic density is due to the heavy Higgs s-channel exchange, allowing for a larger range in \x. 

3.3 Examples of (//, M 2 ) Planes 

We display in Fig. 3 contours of (a,b) the spin-independent and (c,d) the spin-dependent 
elastic scattering cross sections, in the cases tan/5 = 10, and (a,c) m = 100 GeV and 
m A = 500 GeV, (b,d) m = 300 GeV and m A = 500 GeV. Because the cross sections vary 
relatively rapidly, we have not included all the decade cross-section contours in panels (c,d). 

We see again in this case the suppression in the spin-independent cross section for jj, ~ 
— 100 GeV and M 2 ^ 200 GeV, apart from which the cross sections decrease with increasing 
at least within the GUT stability region. The spin-dependent cross section, on the other 
hand, starts rising again at large reflecting the fact that one has traversed a cancellation 
in the scattering matrix element of the type described in the previous subsection. Once 
again, the cancellation is found between the 10~ 7 pb contours just outside the GUT stability 
curve. 

For the values of parameters chosen in Fig. 3 (a,c), we obtain a spin-independent cross 
section a S i = 2 x 10 -9 pb and a spin-dependent cross section a S D = 3 x 10~ 6 pb in the 
CMSSM (cf. the position of the crosses in the figures). In the NUHM, we find that the 
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tan (3 = 10 , m 1/2 = 300 , m = 100 tan (3 = 10, m 1/2 = 300 , m = 300 
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Figure 2: Contours of (a,b) the spin-independent and (c,d) the spin- dependent cross sections 
(solid black lines) are superimposed on the NUHM (fi,m A ) planes for tan/5 = 10, (a,c) 
m = 100 GeV and m 1/2 = 300 GeV, (b,d) m = 300 GeV and m 1/2 = 300 GeV. The 
shadings and line styles are the same as in Fig. 1. Here, we see in addition the constraint 
from b —± 37. The excluded region is medium (green) shaded. The crosses denote the CMSSM 
points for these choices of mo and m^. 
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range of possible cross sections for this case is 2.5 x 10 n pb < asi < 9 x 10 8 pb and 
1.5 x 10" 9 pb < d 5D < 2 x 10~ 3 pb. For the parameters in Fig. 2 (b,d), the CMSSM spin- 
independent cross section is asi = 3 x 10~ 9 pb, while the spin-dependent cross section is 
<y sd = 10~ 5 pb. Note, however, that this CMSSM point would be excluded due to a excessive 
value for VLh 2 ( > 1). In the NUHM, this parameter choice is allowed and gives the range 
10~ 13 pb < a S i < 2 x 10" 8 pb and 9 x 10~ 9 pb < a SD ^ 2 x 10" 4 pb for the elastic cross 
sections. 

4 Allowed Ranges of Elastic Cross Sections 

Following our discussion in the previous Section of some important features in a few exam- 
ples of parameter planes in the NUHM, we now display the ranges of elastic scattering cross 
sections permitted by various theoretical and experimental constraints. We start with the 
specific NUHM parameter planes discussed above, and then go on to generalize the discus- 
sion. In each of the specific planes, we show the effect on the allowed cross section when the 
phenomenological and cosmological constraints are applied successively. We start with the 
very simple requirement that the LSP be a neutralino with m x ± ^ 103.5 GeV and m g ^ 99 
GeV. We then apply either the Higgs cut or the b — > 57 cut. Our standard cut is defined to 
include these two in addition to the appropriate value for flh 2 . Following the standard cut 
we apply sequentially the GUT stability constraint and the constraint due to — 2. 

4.1 Specific Planes 

We consider first the (mi/ 2 ,m ) plane for tan/5 = 10 and /i = m A = 400 GeV that was 
displayed earlier in Fig. l(a,c). The horizontal axes in the various panels of Fig. 4 correspond 
to the LSP mass m x , and the vertical axes show the ranges of (a,c,e) the spin-independent 
and (b,d,f) the spin-dependent elastic cross sections. The first row of panels (a,b) shows 
the ranges allowed by our cuts on the LSP (dark lines), m h (lighter lines) and standard cut 
(shaded), the second row (c,d) displays the further impact of the GUT stability constraint, 
and the third row (e,f) implements all the cuts, including that on — 2. 

We see in Fig. 4(a,b) that the scattering cross sections rise monotonically with m x , 
except close to the upper limit on m x , reflecting the increase with mi/ 2 that was already 
commented in connection with Fig. 1. We also note that the ranges of cross sections allowed 
for any fixed value of m x is very restricted, reflecting the fact the the contours of equal cross 
section in Fig. l(a,c) are almost vertical in the parameter range of interest. The break that 
appears in the middle of Fig. 4(a,b), and is seen more clearly in (c,d), reflects the range of 
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tan p = 10 , m = 100 GeV, m A = 500 GeV tan (5 = 10, m = 300 GeV, m A = 500 GeV 




H (GeV) n (GeV) 



Figure 3: Contours (a,b) of the spin-independent and (c,d) spin- dependent cross sections 
(solid black lines) are superimposed on the NUHM (//, M 2 ) planes for tan/3 = 10, (a,c) 
m = 100 GeV and m A = 500 GeV, (b,d) m = 300 GeV and m A = 500 GeV. The shadings 
and line styles are the same as in Fig. 2, and the crosses denote the CMSSM points for these 
choices of m and rriA ■ 
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Figure 4: Cross sections allowed in a slice through the NUHM parameter space provided by 
the (mi/2,mo) plane for tan/? = 10 and fi = tua = 400 GeV. Panels (a,b) are with our cuts 
on the LSP (dark lines), nih (lighter lines) and standard cut (shaded), panels (c,d) also have 
the GUT stability constraint imposed, and panels (e,f) apply all cuts, including the possible 
g^ — 2 constraint. The right (left) panels show the spin- (in) dependent cross section, plotted 
against m x . 
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mi/2 ~ 500 GeV where Q x h 2 is suppressed below the preferred cosmological range by rapid 
direct-channel XX annihilation via the H, A poles. 

The most relevant effect of the extra GUT and g^ — 2 constraints is to reduce the range 
of mi/2 and hence m x . As we see in Fig. 4(c,d), the GUT stability constraint removes the 
points in this NUHM parameter plane that have the largest elastic scattering cross sections. 
Finally, as we see in Fig. 4(e,f), the — 2 constraint confines our attention to points in the 
NUHM parameter below the rapid-annihilation channel. In this particular case, there is a 
narrow preferred range of the spin-independent cross section around 3 x 10~ 9 pb, and the 
preferred range of the spin-dependent cross section is around 6 x 10~ 6 pb. 

Our second example is the (/IjUIa) plane for tan/3 = 10,mi/ 2 = 300 GeV and m = 
100 GeV, displayed previously in panels (a,c) of Fig. 2. As seen in Fig. 5(a,b), when one 
imposes the standard Q x h 2 , m h and b — > constraints, the cross sections generally decrease 
with m x . However, in the spin-dependent case the cross section reaches a locus of near 
zeroes, after which it rises again. These cancellations are avoided in this case when the 
Higgs and b — > 37 constraints are applied. As then seen in Fig. 5(c,d), the GUT stability 
constraint, which removes portions of this NUHM parameter plane at large and small 
rriAi strengthens in this case the lower bounds on the cross sections. The effect of the g^ — 2 
constraint is less marked in this case, as seen in Fig. 5(e,f). The final allowed ranges of the 
cross section are considerably wider than in the previous example: ~ 10~ 9 to ~ 10~ 8 pb in 
the spin-independent case and ~ 10~ 7 to ~ 10~ 4 pb in the spin-dependent case. 

Our final example is the (//, M 2 ) plane for tan (3 = 10, m = 100 GeV and itla = 500 GeV, 
shown in Fig. 6. In this case, the standard cuts allow a particularly wide range of cross 
sections, varying by infinite (over 3) orders of magnitude for the spin- (in) dependent case. 
The different regions allowed by the standard cuts reflect the different branches of parameter 
space in Fig. 3. We note that some of these branches are due to the differences between 
positive and negative /i. We also note that some of the boundaries are due to our imposed 
cutoff of < 2 TeV. In particular, had we allowed for larger values of we would 
have found larger neutralino masses, and the lower bounds for both the Higgs and b — > S7 
cuts (which differ for positive and negative ji) would also have been lowered. The general 
tendency of the standard cuts is to decrease the cross sections with increasing m x , though 
with considerable variation. The GUT stability constraint in this case removes a region at 
large M 2 and hence m x , that removes the points with the lowest cross sections. This effect 
is particularly marked for the spin-dependent case, where the range is now reduced to (a 
mere) 5 orders of magnitude. The — 2 constraint further raises the lower bounds on the 
cross sections, so that they vary through just over 3 (under 2) orders of magnitude in the 



13 



m m = 300 GeV, m„ = 100 GeV, tan R = 10 



m,„ = 300 GeV, m„ = 100 GeV, tan R - 10 





Figure 5: Cross sections allowed in a slice through the NUHM parameter space provided by 
the (/i,mi) plane for tan/5 = 10 and mu2 = 300 GeV, mo = 100 GeV. The selections of 
points in the different panels are the same as in Fig. 4, as are their indications. In addition, 
panel (a) also shows the b — > s^y cut (dashed lines). 
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spin- (in) dependent case. 

Notice that there is an isolated point at m x ~ 80 GeV. This corresponds to a very narrow 
region just to the right of the chargino mass-bound line in Fig. 3(a,c) at /j, ~ 120 GeV, M 2 ~ 
240 GeV, which is not visible because of the plotting resolution. What happens is that as // 
decreases, m x falls below the threshold, and hence the annihilation cross section decreases 
leading to an acceptable value of Q x h 2 . However, as \x decreases further, neutralino-chargino 
coannihilation becomes stronger, suppressing Q x h 2 again to be less than 0.1. 

This brief survey shows the importance of implementing correctly the GUT stability 
constraint, which may (in different cases) bound the cross sections either above or below. As 
many authors have previously pointed out in the CMSSM case, the — 2 constraint is also 
potentially important. In certain cases, it can also strengthen significantly the lower limits 
on the NUHM cross sections. 

4.2 General Analysis 

Equipped with the above information about some specific examples, we now make a general 
analysis of the possible values of the elastic scattering cross sections. In the first place, we 
concentrate on the case tan/5 = 10, but relaxing the previous restricted choices of other 
parameters that we took as examples. To produce the plots, we generate random points 
(about 30000 points for each plot) within the following ranges: 

100 GeV < m V2 < 1500 GeV, 
< m < 1000 GeV, 
-2000 GeV < /i < 2000 GeV, 

90 GeV <m A < 1500 GeV. (6) 

We first impose the same standard experimental and phenomenological constraints discussed 
earlier, namely: a consistent electroweak vacuum, 0.1 < Vl x h 2 < 0.3, m h > 114 GeV and the 
b — > S7 constraint. In a previous paper on the CMSSM [3], we rescaled the elastic scattering 
cross sections for models that predicted fl x h 2 < 0.1 by the factor il x h 2 /0.1, so as to account 
for the fact that the neutralino could constitute at most this fraction of the galactic halo. 
In the results that follow, we show how this rescaling affects the upper limits on the cross 
sections. Next we apply the GUT stability cut, and finally we show the implications of 
imposing the cut on — 2. The resulting ranges of the elastic scattering cross sections for 
tan (3 = 10 are displayed in Fig. 7. 

We note that the spin-independent cross section shown in Fig. 7(a) may be as large as 
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m A = 500 GeV, m„ = 100 GeV, tan p = 10 




m A = 500 GeV, m = 100 GeV, tan P = 10 



100 200 300 400 500 600 



(GeV) 



m A = 500 GeV, m = 100 GeV, tan P = 10 




100 200 300 400 500 600 



m x (GeV) 



m A = 500 GeV, m = 100 GeV, tan P = 10 




m x (GeV) 




Figure 6: Cross sections allowed in a slice through the NUHM parameter space provided by 
the (/i, M 2 ) plane for tan/3 = 10 and m = 100 GeV, rriA = 500 GeV. The selections of 
points in the different panels are the same as in Fig. 5, as are their indications. 
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100 200 300 400 500 600 



m x (GeV) 

Figure 7: Ranges of (a) the spin-independent and (b) the spin- dependent cross sections for 
tan/? = 10. The ranges allowed by the standard cuts on Q x h 2 , rrih and b — *■ sj have dark 
shading, those still allowed by the GUT stability cut have medium shading, and those still 
allowed after applying all the cuts including — 2 have light shading. The pale shaded region 
corresponds to the extra area of points with low relic densities, whose cross sections have been 
rescaled appropriately. 
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a few xlO -8 pb, decreasing only slightly as m x increases, whilst values lower than 10~ 13 pb 
cannot be excluded. The raggedness of the upper limit on the cross section reflects the fact 
that our sampling produced very few points with such large cross sections: values in between 
the crags cannot be excluded, but must be very rare. Between the crags there are valleys, 
below which the density of points is significantly larger. The lowest values of the cross 
section occur for /i < 0, where cancellations are possible in the spin-independent scattering 
matrix element, as discussed in a previous paper [2]. The GUT stability constraints exclude 
some low cross-section values at both small and large m x , but do not provide an overall 
lower bound. It does however, lower the upper bound by as much as a factor of about 5. 
Negative /x, and hence very low cross-section values, would be excluded by the putative — 2 
constraint, as shown in Fig. 7(a). For \i > 0, we find (not shown) spin-independent cross 
sections only above 1CT 11 pb. 

Fig. 7(a) also displays the region (pale shaded) which survives all standard cuts, except 
that Vth 2 < 0.1. For points in this region, we have rescaled the cross section by a factor of 
fi/i 2 /0.1, to allow for the fact that the LSP could not in this case make up all the cold dark 
matter in the Universe, and hence a fortiori in the galactic halo. As one can see, many of the 
spaces between the crags are now filled in with such points, but very few give significantly 
larger cross sections. 

In the case of the spin-dependent cross section shown in Fig. 7(b), the upper limit is better 
defined, and decreases monotonically from ~ 3 x 1CT 3 pb for m x ~ 80 GeV to ~ 5 x 10~ 6 pb 
for m x ~ 650 GeV. Cross sections lower than ~ 10~ 10 pb are possible for either sign of 
/x, even after imposing the GUT stability cuts. In this case, the points with rescaled cross 
sections enhance the cross section by a factor of about 3 at low neutralino masses. 

The isolated point in both panels at m x ~ 80 GeV now corresponds to a narrow region 
around fx ~ — 110 GeV and large m A > 1000 GeV, which is between the m h = 114 GeV and 
m x ± = 103.5 GeV lines. Its existence is very sensitive to the implementation of the Higgs 
mass bound. 

We note that the ranges allowed by g M — 2 are relatively restricted. For a start, we find 
that 110 GeV ^ m x < 200 GeV. Moreover, even within this range, very low cross-section 
values are excluded. Overall, we find ranges between ~ 2 x 10~ 8 pb and ~2x 10~ 10 pb for 
the spin-independent cross section, and between ~2x 10~ 4 pb and ~2x 10~ 8 pb for the 
spin-independent cross section. 

Ranges for the spin-independent and -dependent cross sections for tan f3 = 20 are shown 
in Fig. 8. Looking first at the spin-independent cross section in panel (a), we see that 
our standard cuts on rrih, Q. x h 2 and b — > S7 would allow somewhat larger values than for 
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tan (3 = 10. This difference is less marked when the GUT stability cut is also applied, except 
for some exceptional parameter choices at small m x . The jaggedness of the peaks is more 
pronounced at this value of tan/5. Once again, we emphasize that while we do not expect 
the area between the peaks to be empty, the density of points there is extremely low. When 
one keeps the low-relic-density points (rescaled appropriately), we see that indeed the crags 
are filled in to some extent. For tan/3 = 20 these points do not enhance the cross section 
significantly. The values of m x allowed by — 2 are larger for tan (3 = 20 than for tan f3 = 10, 
and lower cross sections are also attainable. Overall, the spin-independent cross section may 
vary between ~3x 10~ 7 pb and 10~ 10 pb when tan/3 = 20. 

In the case of the spin-dependent cross section shown in Fig. 8(b), the upper and lower 
bounds for tan (3 = 20 are very similar to those for tan /3 — 10 if only the standard and GUT 
cuts are applied. However, slightly lower values of the cross section are allowed when the 
g^ — 2 constraint is applied, without a strong dependence on the value of m x (which may be 
somewhat larger than in the case of tan/3 = 10). For tan/3 = 20, the spin-dependent cross 
section ranges between ~2x 10~ 3 pb and ~ 10~ 8 pb. 

Continuing now to tan (3 = 35, as shown in Fig. 9, we see again that the spin-independent 
cross section may be somewhat larger still than for tan/3 = 20, though the same is not true 
for the spin-dependent cross sections. Once again, the <? M — 2 constraint allows larger values 
of m x as tan (3 is increased, opening up the possibility of a smaller cross section, particularly 
in the spin-dependent case where a cancellation may occur, potentially suppressing the cross 
section by a couple of orders of magnitude. Both these tendencies are accentuated in the 
case tan/3 = 50, as shown in Fig. 10. For both tan/3 = 35,50, the spin-independent cross 
section may, in some exceptional cases, rise above 10~ 6 pb, even after implementing all the 
cuts. It may also drop as low as 10~ 10 pb. The maximal spin-dependent cross section is 
above 10~ 4 pb in the tan/3 = 35 case, and somewhat below 10~ 4 pb in the tan/3 = 50 case. 
In both cases, it may also drop as low as 10~ 10 pb. 

Finally, in Fig. 11, we display the allowed ranges of (a) the spin-independent and (b) the 
spin-dependent cross sections when we sample randomly tan (3 as well as the other NUHM 
parameters. We see no big surprises compared with the previous plots for individual values 
of tan /3, but observe that the boundaries of the shaded regions are very ragged, reflecting 
the finite sample size. After incorporating all the cuts, including that motivated by — 2, 
we find that the spin-independent cross section has the range 10~ 6 pb ^ a S i ^ 10~ 10 pb, and 
the spin-dependent cross section has the range 10~ 4 pb ^ a S D ~ 10~ 10 pb, with somewhat 
larger (smaller) values being possible in exceptional cases. If the g M — 2 cut is removed, the 
upper limits on the cross sections are unchanged, but much lower values become possible: 
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Figure 8: Ranges of (a) the spin-independent and (b) the spin- dependent cross sections for 
tan/? = 20. The shadings are the same as in Fig. 7. 
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a S i < 10~ 13 pb and a SD < 10~ 10 pb. 

5 Conclusions 

We have discussed in this paper the possible ranges of spin- (in) dependent elastic cross sec- 
tions in the NUHM, in which the input soft supersymmetry-breaking masses of the Higgs 
doublets are allowed to differ from those of the squarks and sleptons - which are still as- 
sumed to be universal. Fig. 11 summarizes the results, including the flexibility of varying 
tan/5 as well as the input scalar and fermion masses m ,i/2, the Higgs mixing parameter \i 
and the the pseudoscalar Higgs mass m A . In this paper, we have not varied the trilinear 
soft supersymmetry-breaking parameter A, whose effective low-energy value has in any case 
a limited range when renormalized from the GUT scale downwards, so that it does not have 
a large effect on the cross sections we study. 

We have stressed in this paper the importance of incorporating consistently all the avail- 
able phenomenological constraints from laboratory experiments and cosmology. We have also 
stressed the importance of taking into account the running of the NUHM parameters over the 
full range of scales between the GUT scale and the electroweak scale. As we have discussed 
explicitly, the requirement that the effective scalar potential be stable at the GUT scale 
restricts the allowed variations in the non-universalities of the soft supersymmetry-breaking 
Higgs masses. 

The effects of our phenomenological cuts and this GUT stability requirement can be 
seen in Fig. 11, and also in the previous Figs. 7, 8, 9 and 10 for tan/3 = 10,20,35 and 50, 
respectively. Some examples of our analysis for specific slices through the NUHM parameter 
space can be seen in earlier figures. 

In general, we find that cross sections may differ by a few orders of magnitude from those 
found in the CMSSM, in which the soft supersymmetry-breaking Higgs masses are assumed 
to be universal with the slepton and squark masses at the GUT scale. However, the spin- 
independent cross section normally lies well below the present experimental sensitivity. Only 
in a few exceptional cases do we find a cross section as large as the present experimental 
sensitivity asi ~ 10~ 6 pb, and a sensitivity asi ~ 10~ 10 pb would be required to cover 
most of the preferred domain of NUHM parameter space. Even this sensitivity would be 
insufficient if one disregards the indication from — 2, which is the only constraint that 
motivates a useful upper bound on the sparticle mass scale, and hence a useful lower bound 
on a si- In the case of the spin-dependent cross section, values of a SD as low as ~ 10~ 10 pb 
cannot be excluded even if one takes seriously the — 2 constraint. 
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The next logical step in the exploration of the MSSM, relaxing further the assumption of 
full scalar-mass universality as in the CMSSM, is to allow the soft supersymmetry-breaking 
slepton and squark masses to differ at the GUT scale. The allowed ranges of the effective 
low-energy slepton and squark masses, after renormalization, will be restricted by analogues 
of the GUT stability constraints we have applied in this paper. In particular, we note that 
general choices of the effective low-energy slepton and squark masses may lead (in particular) 
to tachyonic squarks below the GUT scale, when renormalized to higher scales. We will 
explore in future publications the effects on the parameter space and the ranges of cross 
sections of applying consistently the GUT stability constraints to the general non-universal 
MSSM. 
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